Self-organization and Mechanical Properties of Active Filament Bundles 
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A general framework for the description of active bundles of polar filaments is presented. The 
activity of the bundle results from mobile cross-links, that induce relative displacements between 
the aligned filaments. Our generic description is based on momentum conservation within the 
bundle. By specifying the internal forces, a simple minimal model for the bundle dynamics can be 
derived, capturing a rich variety of dynamic behaviors. In particular, contracted states as well as 
solitary and oscillatory waves appear through dynamic instabilities. We present the full bifurcation 
diagram of this model and study the effects of a dynamic motor distribution on the bundle dynamics. 
Furthermore, we discuss the mechanical properties of the bundle in the presence of externally applied 
forces. Our description is motivated by dynamic phenomena in the cytoskeleton and could apply 
to in vitro experiments as well as to stress-fibers and to self-organization phenomena during cell- 
locomotion. 
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I. INTRODUCTION 

The cytoskeleton of eucaryotic cells is a complex three 
dimensional network of protein filaments, most promi- 
nently actin filaments and microtubules Q, Q . Its elastic 
and viscous properties are essentially defining the me- 
chanical or material properties of living cells. This net- 
work resembles in many aspects a polymer solution or 
a gel. The main difference from usual polymer mate- 
rials is its intrinsic activity. In fact, the cytoskeleton 
is constantly remodeled through the polymerization and 
depolymerization of filaments, as well as through the 
formation and breakup of cross- links. In addition, the 
cross- links may be active, leading to further dynamics. 
Active or mobile cross-links are provided for example by 
molecular motors which are specialized enzymes which 
transduce the chemical energy of a fuel to motion along 
filaments Q, 0, 0, 3 • All these activities are regulated by 
the cell which is thus able to direct intracellular trans- 
port, to separate its chromosomes and to cleave during 
cell division, to exert forces on the environment, or to 
move on a substrate. 

The study of active polymer systems requires com- 
pletely new tools and techniques as compared to the well 
developed analysis of equilibrium properties, that relies 
on powerful concepts of equilibrium statistical physics. 
Indeed, such systems are intrinsically far from equilib- 
rium, and the dynamics at equilibrium, that is usually 
studied in polymer physics [5|, |fj, ML is not sufficient for 
the description of active systems 8]. On the contrary, 
experimental studies of the cytoskeleton under simpli- 
fied conditions have revealed its ability to self-organize. 
Namely, the contraction of filament bundles |jj , the for- 
mation of asters and vortices 0, 0, 0, ^| , as well as 
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the formation of networks [bH ] were found in vitro. Using 
a cell extract, even the formation of bipolar spindles with- 
out microtubule organizing centers has been seen |l5j| . 
Cell fragments containing only the actin cytoskeleton, 
but neither the nucleus nor microtubules, can propagate 
on a substrate |l6(, where the locomoting state coexists 
with a stationary spherically symmetric state 01 • I n 
a mixture of actin filaments and myosin molecular mo- 
tors, active reptation in a polymer solution has been ob- 
served |l8j . Let us finally mention, that experiments 
probing mechanical properties of living cells have re- 
vealed active responses of the cytoskeleton to external 
forces, see, e.g., Ref. |l9| . 

First steps towards a theoretical understanding of ac- 
tive polymer systems have mostly aimed at describing 

pattern-formation. In one-dimensional filament bundles, 
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polarity sorting |2f|, contraction 21, 22], and propagat- 
ing waves (2^1 have been observed. Self-organization has 
also been seen to induce bending waves and complex mo- 
tion in axonemes |2J, |25| . In higher dimensions, the ef- 
fects of active cross-links on the formation of orientation 
patterns in svstems of spatially fixed filaments have been 
studied |2til27| and the generation of filament currents by 
active cross-links has been discussed .2£|. Furthermore, 
the viscoelastic response of solutions of semi- flexible poly- 
mers and active centers has been studied pfij . 

Active filament bundles provide very simple examples 
of active filament networks and can be discussed by a 
one-dimensional description. Note, however, that in ad- 
dition to their simplicity, such filament bundles actually 
occur in animal cells. They are, for example, part of 
stress fibers which generate contractile forces, and of the 
contractile ring in dividing cells [2,]. In vitro, the contrac- 
tion of actin bundles in the presence of myosin motors has 
been observed 0. 

Motivated by the dynamics of the cytoskeleton, we de- 
velop a general framework based on momentum conser- 
vation, to describe the physics of bundles of aligned fil- 
aments in the presence of active cross-links. Since actin 
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filaments and microtubules are polar as they have two 
structurally different ends, we consider polar filaments. 
The cross-links are mobile and considered to be formed 
by small aggregates of molecular motors of one type. We 
discuss simple scenarios in order to study dynamic phe- 
nomena and mechanical properties of such systems. A 
minimal model which has been introduced in earlier pub- 
lications |22l |23| can be derived in our general framework 
using approximations and simplifications. This model al- 
ready exhibits many of the phenomena that occur in such 
systems. However, it neglects changes in the distribu- 
tion of motors due to the dynamics of the system. This 
dynamics of the motor distribution can be taken into 
account within our general framework. Furthermore, we 
discuss the mechanical properties of active bundles in the 
presence and absence of external forces. 

The outline of our manuscript is as follows. In sec- 
tion [n] we introduce the general description for active 
filament bundles based on momentum balance. Using 
this formalism we derive dynamic equations for the sys- 
tem. In Section IIIII we show how the minimal model 
can be derived from the general equations. We review its 
properties and extend previous work towards a discussion 
of the full bifurcation diagram. In Section ITvl we study 
the effects of the dynamics of the motor distribution on 
the filament dynamics. The active mechanical properties 
of a filament bundle can be derived systematically in the 
framework introduced in section [H] We discuss the bun- 
dle mechanics in Section [V] and study tense states bal- 
anced by external forces applied at the ends. The paper 
concludes with a discussion of our results in Section IVT1 
which relates our theoretical framework to experimental 
situations. The appendices contains a detailed analysis 
of bifurcations in the minimal model as well as the effects 
of filament adhesion to a substrate. 



II. PHENOMENOLOGICAL DESCRIPTION OF 
ACTIVE FILAMENT BUNDLES 

We introduce a general description for the dynamics of 
active filament bundles. The bundle is described in one 
dimension using density profiles of filaments and motors. 
The dynamics of these densities is governed by currents 
which are generated by interactions between filaments 
and motors. Dynamic equations can be derived on the 
basis of momentum balance. This general procedure can 
be carried out most conveniently using simplifying as- 
sumptions. In particular, we assume a low motor den- 
sity or low duty ratio of motors such that interactions 
between filament pairs dominate, a local friction of fila- 
ments with the environment and we neglect the possibil- 
ity of passive cross-linkers. Furthermore, we assume that 
filament lengths remain fixed, i.e. we neglect polymeriza- 
tion and dcpolymerization of filaments, and we assume 
that filaments cannot change their orientation. However, 
many of the qualitative behaviors displayed by the re- 
sulting equations are more general and are also found 



in purely phenomenological descriptions which are not 
based on these simplifying assumptions [3(|. The dy- 
namic equations we discuss here represent a mean field 
theory of filament bundles where fluctuations do not ap- 
pear explicitly but give rise to diffusive terms. 

A. Densities of filaments and motors 

The bundle is characterized by the number densities 
of filaments and of motor complexes projected on the 
bundle axis, which leads to an effective one-dimensional 
description. Since filament bending and entanglements 
can be ignored in the bundle, we describe filaments as 
rigid rods. Filaments are aligned along the x-axis and 
we distinguish the two sub-populations of filaments with 
their plus-end pointing into the positive and negative x- 
direction, respectively. These populations are described 
by the densities c + and c~, such that, e.g., c + (x)dx gives 
the number of filaments with their plus-end in the pos- 
itive x-direction and their center located in the interval 
[x, x + dx]. We assume, that motors are small compared 
to the filament length and will be treated as point-like 
in our description. The number density of motors is de- 
noted by m. 

The filament and motor densities satisfy the conserva- 
tion laws 

d t c+ = Dd 2 x c+-d x j+ (i) 

d t c~ = Dd 2 x c~ - d x J~ (2) 
d t m = D m d x m — d x J (3) 

Here, the currents J ± and J are generated by the active 
interaction between motors and filaments. The densities 
c + and c~ are conserved separately since we do not al- 
low filaments to change their orientation. Fluctuations 
in the system give rise to diffusive terms with diffusion 
coefficients D and D m of filaments and motors. While 
the diffusion of motors could be expected to result from 
thermal fluctuations, the diffusion of filaments is gen- 
erated effectively by fluctuations of the forces induced 
by motor-filament interactions. For long filaments, the 
contribution of thermal fluctuations to the diffusion co- 
efficient D is negligible. We return to this point in Sect. 

eh 

B. Momentum balance 

In the absence of external forces, the total momen- 
tum is conserved in the filament bundle. Forces acting 
within the bundle lead to an exchange of momentum with 
the environment or between filaments. Since filaments 
are treated as rigid, extended objects with momentum 
distributed along the full length of the filament, we in- 
troduce the momentum densities tt (x,y). These densi- 
ties represent the momentum at position y carried by all 
plus- or minus-filaments, respectively, with their centers 
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located at position x. The momentum balance can then 
I 

d t n+(x,y) + d y (T + (x,y) - f£ t (x,y) 
d t ir~(x,y) + d y (T~{x,y) - f-~ t {x,y) 



Here, momentum flux along filaments centered at x is 
given by the tensions a (x,y). Momentum exchange 
between filaments is nonlocal and described by the in- 
ternal force densities f^ t (x,y), which include all active 
filament interactions via motors. The force densities 

fa( x ,y), fm( x ,y)> and fia( x ,y) are source and sink 
terms, describing momentum exchange with the environ- 
ment. They result from friction with the fluid (fa), from 
motors moving along a single filament (/ m ), and from ex- 
ternal forces (/cxt)- Here, x refers to filaments with center 
at position x, while y denotes a position in space, where a 
force is acting and momentum is exchanged. Momentum 

I 



Incrtial terms are negligible in a slowly moving bundle, 
such that we can set dtir^ — 0. Equations and 
then express a balance of forces. 

In the most simple case where friction is local, we can 
write for the density of friction forces 

f±(x,y) = V J ± (x)R(x-y) (8) 

Here, 77 is a friction coefficient per unit length and R(x) is 
a function characterizing the distribution of energy dissi- 
pation along moving filaments. If all filaments are of the 
same length £, a simple choice is R(x) = 1 for \x\ < £/2 
and R(x) — otherwise. However, the function R(x) can 
also account for situations with a distribution of filament 
lengths. Then, R(x) is related to the probability, that a 
given filament is longer than 2\x\. 

The forces exerted by motors moving along a single 
filament are linear in the filament and motor densities, 

fm ( x , V) = T77mrm(y)c ± (x)R(x - y) , (9) 

where 7y m is the friction coefficient corresponding to single 
motors and the coefficient T characterizes the binding to 
and motion on filaments of individual motors 



be expressed as 



= ff[(x,y) + f+(x, y ) + f+ t (x,y) (4) 

= fa( x -,y) + fm( x ,y) + /ext(^y) ( 5 ) 
1 

conservation in the absence of external forces requires, 
that 

J dx[f+ t (x,y) + fr t (x,y)]=0 . (6) 

This implies, that any force generated by active cross- 
links on a filament is balanced by an opposite force act- 
ing on other filaments, i.e., internal forces at a point y 
are balanced when integrated over all filaments. There- 
fore, total momentum H = J dx dy(n + + 7r~) changes 
according to (ignoring boundary terms): 



(7) 

I 

C. Currents of filaments and motors 

While the internal forces are balanced at a point 
y when integrated over all filaments, the total force 
J dyf int (x, v) ac f m g 011 filaments centered at a given posi- 
tion x does not vanish in general. Integration of Eqs. (@J 
and JSJ with respect to y reveals that this force is bal- 
anced by friction forces: 

^{x) = ~ J dy y) + f£(x, y)] (10) 

where £ = J dxR(x) is the average filament length and 
where we have assumed / ox t = 0. Since friction of mo- 
tors is small as compared to filament friction, rj m <c r}£, 
the contribution of can be neglected in most practi- 
cal cases, in particular, for low motor densities. In the 
following, we therefore set ?/ m — 0. 

In order to write explicit expressions for the currents, 
we need a model for the internal forces in the bundle. We 
consider the case when clusters of three or more cross- 
linked filaments form rarely enough, such that their con- 
tribution to the internal forces can be neglected. This 
holds in the case of a low motor density or for motors 
with a low duty ratio, which is the fraction of time a 
motor spends attached to a filament [3l| • If interactions 
between filament pairs dominate, we can split the inter- 
nal forces into those between filaments of the same and 



dt 



n= / dx dy [f+t(x,y) + f eXb (x,y) + f£{x,y) + / fl (x,y) + f+(x,y) + f m (x,y)] 
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those of opposite orientation. We write 

/int = /int + + /int (H) 

and analogously for f~ t . A motor may link two filaments 



and thus exert forces of opposite sign on each of them, 
whenever they overlap. Assuming that the probability 
for two filaments to interact increases quadratically with 
filament density, we write 



J 



fitt ( x > v) = I dz c± ( x ) c± ( z ) R (y ~ x ) R (y - z ) m (y) F±± (z-x,y-x) 



(12) 



and corresponding expressions for f^ 1 . Here, () 
is the average force acting on plus-filaments at a distance 
£ from the center exerted by motors, that interact with 
other minus-filaments located at a distance £ from the 
first, see Fig. ^ The essential feature of motor- filament 
interactions is that the direction of the force applied by a 
motor on a filament is uniquely determined by the orien- 
tation of the filament . The product R(y — x)R(y — z) 
gives the probability, that a filament at x has an over- 
lap at y with a filament at z. Here, the position of a 
filament is given by the position of the filament's cen- 
ter. Analogous expressions hold for the internal forces 
between minus-filaments as well as between filaments of 
opposite orientation. 

The forces F ±± and F ±=F obey the following symmetry 
relations. Momentum balance demands that under an 
exchange of filaments, the force changes sign: 



F±± (Z,0 



-F (- 








(13) 
(14) 



see Fig. Using relation (|12|l , the internal forces satis- 
fying the above equations verify Eq. ©, which assures 
momentum conservation. Space inversion symmetry re- 



quires 



-F~ 
-F~ 



-0 
-0 



(15) 
(16) 



Momentum conservation also determines the non- 
diffusive motor current J with 



VmJ(y) = / dx[f+(x,y) + f m (x,y)} 



(17) 



Here the forces of individually bound motors are defined 
in Eq. ©. 

The continuity Eqs. CEJ-© for the densities together with 
the defining Eqs. (|1(J|> and (|12fl of filament currents and 
Eq. (|17f) for the motor currents provide the full dynamic 
equations of active filament bundles. The functions F ±=t 
depend on the details of the motor-filament interactions 
and could be modified by further proteins bound to the 
filaments. However, the large scale behaviors of the sys- 
tem do not depend on the detailed form of these func- 
tions. In the following sections, we will therefore make 
simple choices which obey the symmetry relations dis- 
cussed above. 
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FIG. 1: Schematic representation of the forces exerted by 
an active cross-link on two filaments of opposite orientation. 
The centers of the filaments are indicated by the dashed lines, 
whereas the dotted-dashed line marks the position of the ac- 
tive cross-link. The arrows indicate the direction of the forces 
applied by motors on the filaments. 



III. THE MINIMAL MODEL 

The minimal model has been introduced in Ref. [2^ 
as a simple model for filament dynamics. It can be ob- 
tained from the general equations derived in the previous 
section by choosing R(x) = 1 for \x\ < £/2 and R(x) = 
otherwise. The forces are chosen to behave as 



sgn(£) 



(18) 



where sgn(£) = 1 for £ > and —1 otherwise, which 
represents the simplest choice compatible with the sym- 
metry requirements. Furthermore, we assume in the min- 
imal model that the motor distribution is homogeneous 
and its dynamics can be neglected. 

The resulting dynamical equations are most conve- 
niently expressed in dimensionless form. We define x = 
xjl and measure lengths in units of the filament length 
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£, and a dimensionless time variable t = tD/i 2 . Further- 
more, we introduce dimensionless densities c = cl. Sup- 



J 



pressing the tildes, the dynamic equations of the minimal 
model can be written as 



d t c + (x) = d 2 c+(x) -ad x d£[c + {x + t)-c + (x~£)}c + (x) 



+(3d x / d^c-(x + Oc + {x) 



d t c (x) — <9 2 c (x) -ad x j d£ [c (x + f) - c (x — £)]c (x) 



-(3d x j d£c+(x + Oc-(x) 



(19) 



(20) 



where a and (3 are dimensionless coupling constants char- 
acterizing the strength of the motor forces defined in 
Eqs. I]18p. It follows from the dynamical equations, that 
the homogeneous state c^^x) — = const, is a station- 
ary solution for all values of the parameters. 



A. Oriented bundles - Contraction 

If all filaments are of the same orientation, one is left 
with a single equation 

d t c(x) = d 2 x c{x) - ad x f d£ [c(x + 0- c(x - £)] c(x) . 
Jo 

(21) 

Here, c represents either c + or c , depending on the 
orientation of the filaments. This nonlinear integro- 
differential equation is the most simple description of the 
active dynamics of a filament bundle. Many of the ba- 
sic physical principles underlying self-organization of fil- 
ament bundles can already be discussed using this equa- 
tion. 



(|2*TJl reads 



d_ 

It 



—c k = - (k 2 - 2ac a (1 - 
= A(fc)c fe , 



cos 



fe)) 



Ck 



(23) 
(24) 



for all k. This relation implies that for aco < k 2 /2(l — 
cos k) the mode Ck decays in time, because then X(k) < 0. 
It follows that the most unstable mode is the one cor- 
responding to the smallest non-zero wave number k — 
2irn/L with n = 1. This can be demonstrated using 
(2£) 2 /2(l - cos ^) < (^f.) 2 /2(l - cos 2p) for all n > 1 
which can be verified by induction using the equivalent 
condition n 2 — n 2 cos(27r/L) — l+cos(27rn/L) > 0. There- 
fore, the homogeneous state is linearly stable as long as 
a < a Cl where the critical value a c is determined by 
X(2tt/L) = 0. Explicitly, 



2tt 2 



c L 2 (l-cos(27r/£)) 



(25) 



The critical value a c is positive and decreases with 
increasing Co and L (for L > 1). Note, that for bundle 
sizes L > 1, we have < a c < oo. 



1. Linear stability 

We consider a system of length L with periodic bound- 
ary conditions and study the stability of the homoge- 
neous state with respect to small perturbations. Periodic 
boundary conditions imply that the bundle forms a ring. 
Such rings appear, e.g., in eucaryotic cells, in the late 
stages of cell division. We represent the filament density 
by a Fourier expansion 

c{x) =5> fc e ifa , (22) 

k 

with k = 2irn/L, n = 0, ±1, . . . and where c_fe = c* k . Up 
to first order in the Fourier components Cfc , the dynamics 



2. Contracted states 

If the homogeneous state is unstable, the system 
evolves to an inhomogeneous steady state. We can calcu- 
late this state by numerically solving the dynamic equa- 
tions or, in the vicinity of the bifurcation, using a sys- 
tematic expansion in Fourier modes. To third order in c\ 
the equation for the steady state d t c = reads 

F{a)c l -G{a)\c l \ 2 c l = Q , (26) 

with F(a) = X(2tt/L) and G{a) given by Eqs. (|B6fl and 
(|B7|) . see Appendix |B] Note, that F(a c ) = 0. Expanding 
F and G at a = a c , we find expressions for the Fourier 
amplitudes c\ and C2, given by Eqs. (|B8|I and (|B9(I . This 
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FIG. 2: Schematic representation of the amplitude |ci| of the 
first spatial Fourier-component of stationary states as a func- 
tion of the coupling strength between equally oriented fila- 
ments a. Presented are the cases of a super- and of a sub- 
critical bifurcation. Solid lines represent stable, dashed lines 
unstable solutions. In both cases, the homogeneous state is 
stable for a < a c and unstable otherwise. In the supercritical 
case, i.e., if F/G > for a > a c , the bifurcating solution ex- 
ists for a > ct c and is stable, while in the other case it exists 
for a < ct c and is unstable. In the latter case, one usually 
finds inhomogeneous solutions coexisting with the homoge- 
neous state in an interval [aid, a c]- 



FIG. 3: The amplitude of the first Fourier-component of sta- 
ble stationary solutions of the minimal model for an oriented 
bundle as a function of the interaction strength a. The aver- 
age filament concentration is Co = 0.7 and system size L — 10. 
The inset presents the non-homogeneous stationary solution 
for a = 1.5 . The scenario shown corresponds to a subcritical 
bifurcation, see Fig. 

For large a we find numerically, that transient states 
consisting of several contracted packets can occur, which 
however decay for long times to a steady state with one 
maximum. Their lifetime increases with increasing a. 



solution represents a localized distribution of filaments, 
i.e., a contracted bundle. 

It follows from Eq. I)26[l. that this contracted steady 
state exists if F{a)JG{a) > 0. Depending on whether 
the ratio F/G is positive for a > a c or for a < a c , the 
bifurcation is supercritical and subcritical, respectively, 
see Fig. From Eq. IjBSjl one deduces that the bifur- 
cation is supercritical for system sizes falling within par- 
ticular intervals for which 4/ (4n — 1) < L < 4/ (An — 3), 

71 = 1,2,.... 

Figure |31 presents numerical solutions of the dynami- 
cal equations using an Euler-algorithm with spatial dis- 
cretization A = 0.1. Displayed is the modulus of the first 
Fourier-component |ci| of the attractors as a function of 
a for a system of length L = 10. A region of coexistence 
extends from a = ad up to a c , where the homogeneous 
state becomes unstable. 



3. Contraction dynamics 

The contraction of the bundle is most conveniently dis- 
cussed in an infinite system using the variance of the fil- 
ament distribution 

1 f°° 

o 2 = - dx x 2 c(x) , (27) 

^ * J — oo 

where N = J_ dx c(x) is the total number of filaments, 
as a measure of bundle contraction. Since the center of 
mass of the distribution is immobile due to momentum 
conservation, we have chosen without loss of generality 
(x) = J_ dx xc(x) — 0. The variance changes in time 
as 



dt 



1 f°° 

— dx x 2 d t c(x, t) 

J — oo 

2 f°° , , v 2 
dx c(x, t) + — 



dx xc(x + £,)c(x) 



dx xc(x — £)c(x) 



2 f°° 
2 - —a d££ dx c(x + £)c(x) 

A Jo J-oo 



(28) 
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The final expression reveals two opposing effects. The 
positive constant describes the spreading of the bundle 
due to fluctuations while the second term takes into ac- 
count the effect of the active interactions. The inter- 
action between parallel filaments tends to contract the 
bundle. Note, that distributions for which da 2 /dt = 
do not necessarily correspond to stationary solutions of 
the dynamics (|2T]l . 



B. Bundles of mixed orientation - Solitary waves 

1. Linear stability 

The linearization of Eqs. (|19f) and l|2UI) around the ho- 
mogeneous state c ± (x) = c = const reads in Fourier- 
representation 



We find that at the bifurcation a solitary wave of the 



d 
dt 



A++ A^ 
A-+ A" 



(29) 



where the elements of the matrix A(k) are given by 



A ±:t (fc) = -k 2 - 2a(cos(fc) - l)cj ± 2i/3k4 
A ±=F (fc) = ±2i/3sin(fc)c t . 



(30) 



For a system of length L with periodic boundary condi- 
tions, the wave numbers are k = 2nn/L with n = 0, 1, . . .. 
In presence of the coupling between the filaments of op- 
posite orientation, the matrix A is not diagonal. The 
stability of the modes is determined by the larger of the 
real parts X(k) of the complex eigenvalues of this matrix 
which are given in App. [XI 

We find again that the mode with the smallest wave 
number k = 2t:/L is most unstable and there exists a 
critical value a c where the homogeneous state becomes 
linearly unstable. Furthermore, a c > 0, independently 
of the values of the other parameters, see App. [X] The 
critical value 



c V c 



(31) 



form c ± (x,t) 



: (x — vt) occurs. From momentum 



conservation, it follows that the total filament current 



/ = J dx(J+ + J~) 



(32) 



vanishes. This implies that the propagating filament pat- 
tern is not accompanied by a net filament transport. 
However, as soon as filament adhesion to a substrate is 
introduced, the total filament current associated with a 
solitary wave no longer vanishes and self-organized fila- 
ment transport occurs j23|, see App.IUI 

For weak interactions between filaments of opposite 
orientation, |/3| <C 1, solitary waves can be understood 
intuitively. They emerge from the interaction of a con- 
tracted distribution of filaments of one orientation with a 
homogeneous distribution of filaments of opposite orien- 
tation. This picture suggests a systematic procedure for 
determining solitary waves. Writing c^(x,t) = u ± (a; — 
vt), we can expand u^(x) in powers of (3. For j3 = 0, 
we start from steady states as discussed above, which we 
denote as c ± {x,t) = Uq(x). Solitary waves are obtained 
by assuming that, e.g., Uq is a contracted steady state, 
while Uq is homogeneous. We can now write 



(x) = uf (x) + uf (x)(3 + U2P 2 + 
v = vi{3 + v 3 (3 3 + ... 



(33) 
(34) 



The even terms in the expansion for v vanish by symme- 
try. To lowest order, we obtain 



Vi 



2Cn 



2i sink c Uq k 
— 2a(l — cos k)cQ 



(35) 
(36) 

forfc^O. (37) 



This result agrees with numerical solutions of Eqs. I|19|) 
and (|20|) , see Ref. [2^ • For any /3, solitary wave solutions 
can be obtained near the instability via a systematic ex- 
pansion in Fourier modes, see App. IB 21 



is a function of the remaining parameters, where c = Cg + 
Cq and 8c — Cq — Cq . Here, g is a dimensionless scaling 
function. In some limiting cases, explicit expressions for 
a c can be obtained. For example, for 5c = one finds 
g = Att 2 /L 2 [1 — cos(27r/L)] and in the limit L — > oo 
g = 1 if /3 ^ 0, whereas g = 2/(1 + \Sc\/c) for (3 = 0, 
see Eq. (12511 . These expressions reflect some interesting 
properties of a c , e.g ., it decreases monotonically with 
\Sc\/c, L, and \(3\ [22|. 



2. Solitary waves 

For (3^0, the eigenvalues of A(k) are complex and 
the homogeneous state loses stability through a Hopf- 
bifurcation, leading to solutions that oscillate in time. 



3. Bifurcation diagrams 

The arguments of the previous subsection suggest 
a general scenario of bifurcations as the interaction 
strength a is increased. This scenario emerges from the 
situation with [3 = 0, where the two subpopulations 
evolve independently, if (3 becomes nonzero. For small 
a the homogeneous state is stable. At the critical value 
a c , corresponding to the instability of one filament pop- 
ulation at (3 = 0, a solitary wave occurs. Furthermore, a 
second bifurcation where the solitary wave loses its sta- 
bility occurs at a' c . This second bifurcation is related to 
the point where the second filament population becomes 
unstable for (3 = 0. This bifurcation leads to an oscil- 
lating wave solution which consists of two oscillating dis- 
tributions, which periodically penetrate each other. An 
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FIG. 4: Total filament concentration c = c + +c~ as a function 
of position x and time t for An oscillating wave solution of the 
minimal model for coupling parameters a — 2.5 and (3 — 1., 



average filament densities cj 
size L = 10. 



0.7 and c = 0.3, and system 



example of such an oscillating wave is shown in Fig. 0] 
It is characterized by propagating filament profiles which 
oscillate. 

These arguments allow us to derive, the full bifurcation 
scenario for small values of (3. For ad < a < a c homoge- 
neous distributions and solitary waves coexist. Similarly, 
we find a coexistence of oscillating and solitary waves 
for a' d < a < a' c . Depending on the ratio of plus- and 
minus-filament numbers, Cq /cq, we find different bifur- 
cation scenarios, see Fig. For large ratios, a c < a' d 
and the coexistence regions are separated. In the sec- 
ond case of similar filament numbers shown in Fig. 
a' d < a c and multiple coexistence occurs, in particular 
a coexistence of two different solitary waves indicated as 
Si and 5*2. These waves move into opposite directions. 
For the special case Cg = Cq , the system is symmetric 
under x — > — x and the coexisting waves Si and S% occur 
via spontaneous symmetry breaking. 

For three different values of (3, numerically obtained 
bifurcation diagrams are displayed in Fig. [SJ Shown is 
(l c i~l + l c rl)' i- e -> the time-averaged sum of the ampli- 
tudes of the first spatial Fourier components of both dis- 
tributions as a function of a. For (3 = 0.001 the diagram 
follows closely the curves shown in the bottom panel cor- 
responding to the limit of small (3. As {3 is increased, the 
coexistence regions shrink and the solitary wave £2 disap- 
pears. The bifurcation from solitary waves to oscillating 
waves becomes supercritical for large (3. 
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FIG. 5: Asymptotic solutions of the minimal model rep- 
resented by the time-averaged total amplitude \cf\ + \a[\ 
of the first spatial Fourier-components as a function of the 
interaction strengtha for different values of the interaction 
strength fl. The homogenous state is indicated by H, soli- 
tary waves by S, and oscillatory waves by O. (a)-(f) have 
been obtained by numerical integration of the minimal model, 
where /3 = 0.1 for (a) and (b), /3 = 0.01 for (c) and (d), and 
/3 = 0.001/or(e)and(/). (a), (c), (e) are for c+ = 0.7 and 
Cq = 0.3, while (b), (d), (f) are for Cq = 0.55 and Cq = 0.45 . 
In all cases L — 10. (g) and (h) give a schematic representa- 
tion valid for \/3\ <^ 1 derived from the bifurcation diagram for 
f3 = as explained in the text. Dashed lines indicate unstable 
solutions. 



that the motor distribution remains homogeneous, which 
implies that motors diffuse infinitely fast. We will now 
discuss the effects of the dynamics of the motor distribu- 
tion using the same choices for the motor forces as given 
by Eqs. (|18fl and we assume that filaments are of the 
same length I. 



A. Dynamical equations 



IV. DYNAMIC MOTOR DISTRIBUTIONS 

Motors are actively transported along filaments which 
in general leads to dynamic changes in the motor distri- 
bution |27ll33| . In the previous section, we have assumed 



The dynamic equations, including the dynamics of the 
motor density, have been derived in section II. As in the 
last section, we use dimensionless space and time coor- 
dinates, x — xjl and t — tD/£ 2 as well as dimensionless 
densities c . We introduce a dimensionless motor den- 
sity fh = m/mo, where mo — x Jo ^ x m ( x ) an( i is the 
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system size. Furthermore, we define the dimensionless 
parameters T = T£/D, and D m — D m /D. Suppressing 
the tildes we obtain 

ftc* = dlc^ - - d x J±* (38) 

(39) 



9(771 = D lrL d x m — d x J 



with 



J ±± (x) = A J d£sgn(0 M(x, ^(x + £)c t (x) (40) 
J±T( X ) = tB J d^M{x,S,)c^{x + (,)c ± {x) (41) 

1 

J{x) = r j d£ [c + (x + -c-(x + £)] m(x) (42) 
Here M is the number of motors present in the overlap 



region of two filaments with M(x,Q = f _i, f d( tti(C) 

for £ > and M(x,f) = f*^"? +f dC n»(C) for £ < 0. 
Furthermore, the dimensionless coupling constants A and 
B are related to the parameters a and /? of the minimal 
model and describe effective interactions of a motor with 
a filament pair. 



B. Oriented bundles 

We consider again a system of length L with peri- 
odic boundary conditions and all filaments of the same 
orientation characterized by the density c{x). The dy- 
namic equations linearized at the homogeneous state 
with c(x) — cq and m(x) = mo = 1 read in Fourier 
representation 



d 
dt 



Ck 

m k 



-k 2 + 2A(1 - ^)c Q -2A(cos§ - |sin|)c 2 



-2?Tsin; 



-D m k — iTkco 



m k 



(43) 



where k = 2irn/L, with n being integer. For T/D m —* 0, 
we recover the minimal model discussed in the previous 
section. In contrast to the minimal model, the eigenval- 
ues of (|43|l are complex if T ^ 0. The homogeneous 
state thus loses stability through a Hopf-bifurcation, 
where the most unstable mode occurs at a characteris- 
tic wave- number ko- As a consequence, for larger system 
sizes, solitary waves with multiple maxima appear, see 
Fig. IIV Bl Furthermore, oscillatory waves which can co- 
exist with solitary waves have been found, in contrast to 
oriented bundles in the minimal model. 




2 4 6 8 x 



C. Bundles of mixed orientation 

In the case of mixed bundles, the linearized equations 
read in Fourier representation 



dt ° k 
m k 



A++ A+- A+ m 
A-+ A— A- m 

AUl-f- Am- ^mm 




(44) 



m k 



where the elements of the matrix A(fc) are 

sin k 



A ±± (fc) = 
A ±T (fc) = ±2iB 



k z + 2A{l 
1 — cos k 



—)4±iBk4 



A ±m (fc) 



-2A(cos 



A m± (fc) = T22Tsin 



2 k 
k 



sin - c„ ± 2iBsm — cX c, 



o <-o 



A mm (fc) = -D m k 2 - iTk{c+ - Cq 



FIG. 6: Examples of asymptotic solutions for an oriented bun- 
dle, where the dynamics of the motor distribution has been 
taken into account, (a) shows a solitary wave with a spa- 
tial period of half the system size, (b) an oscillatory wave. 
The solid lines represent the filament distributions, the dashed 
lines the motor distributions, (b) distributions at four differ- 
ent times are shown. They have been shifted in the y-direction 
for better visibility and in the x-direction such that the po- 
sitions of the maxima coincide. The two solutions in (a) and 
(b) exist for the same parameters, which are A — 1.5, T = 1.5, 
Dm = 1., cq = 1., mo = 0.5, and L = 10. 
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FIG. 8: Inhomogeneous steady state of a symmetric system 
with Cq — Cg where the dynamics of the motor distribution 
m has been taken into account. The filament distributions 
are represented by dotted and dashed lines for c + and c~, 
respectively, the full line represents the motor density m. The 
parameters are Cq — Cq = 1, B = T = D n 



1, and L = 10. 



FIG. 7: Regions of stability of the homogeneous state for a 
filament-bundle with dynamic motor distribution. A and B 
characterize the coupling of a filament pair via a motor, where 
the homogeneous state is stable below the lines shown, (a) is 
for the symmetric case Cg = Cg = 1, (b) for Cq = 1 and = 
0.5 . In both cases, L = 10 and D m = 1. Solid lines indicate 
oscillatory instabilities, dashed lines static instabilities. The 
inset presents a case where increasing A can re-stabilize the 
homogeneous state. 



Fig. [7| displays the line in the (A, i?)-plane limiting the 
region of stability of the homogeneous state for different 
values of T. In contrast to the minimal model, it de- 
pends non-monotonically on B. This implies that upon 
increasing B the homogeneous state can be re-stabilized. 
For large enough T we even find reentrant behavior with 
respect to A, i.e., increasing A can re-stabilize the homo- 
geneous state, see Fig. \7\ Furthermore, the figure shows 
that the critical values of A can become negative. This 
implies that the homogeneous state can become unstable 
even if interactions between parallel filaments are absent 
i.e. A = 0. Another interesting observation is that for 
c o ~ c o the homogeneous state becomes unstable with 
respect to stationary inhomogeneous states in certain pa- 
rameter ranges of B. These inhomogeneous stationary 
states represent a new type of solution as compared to 
the minimal model of bundles of mixed orientation. 

Dynamic solutions to this model can be studied nu- 
merically. We find the three types of solutions discussed 
before namely homogeneous states, solitary waves, and 
oscillatory waves. In addition we obtain the inhomoge- 
neous stationary state mentioned above which is a one- 
dimensional analog of asters in higher dimensions and 
consist of localized distributions of motors and filaments, 



see Fig. |HJ In contrast to the minimal model, solitary 
waves can exist even in the absence of interactions be- 
tween filaments of opposite orientation (B = 0) and sta- 
tionary inhomogeneous states can exist for B =/= 0. Inter- 
estingly, the direction of motion of solitary waves can be 
reverted by changing T. 

In summary, the important effects of a dynamic motor 
distribution are (i) the appearance of a characteristic 
wave-length independent of the system size, (ii) the ap- 
pearance of aster-like solutions, and (iii) that the exis- 
tence of inhomogeneous states requires a symmetric sit- 
uation with ct = CrT . 



V. CONTRACTILE TENSION AND EXTERNAL 
FORCES 

So far, we have been focusing on dynamical properties 
of active bundles. We now discuss mechanical proper- 
ties of the bundle such as bundle tension and the role 
of applied external forces. The total tension at a 

point y within the bundle is obtained by integrating the 
contributions of all filaments, i.e., 

E(y) = f dx[a + (x,y)+a-(x,y)] , (45) 

where a (x, y) has been introduced in Eqs. and JSJl. 
Taking into account momentum conservation, Eq. J(|J|, 
one finds in the absence of external forces 

^E(y) =Jdx [/+ + / fl " + /+ + /"] . (46) 

This equation allows us to calculate the tension profile 
by means of Eq. JSJ if the currents J are known. 
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For simplicity, we neglect again f m . The total tension 
in the bundle can be written as 

E(y) =S^(y) + S^(y) . (47) 

I 

E=»(tf) = 

^(y) = 

where -£Q{x) = R(x) with Q(0) = 0, and E^ and E^ J 
are constants of integration. These expressions seem to 
imply that the tension at position y depends on the global 
state of the bundle since Q(x) ^ for all x ^ 0. However, 

I 



and 



In contrast to the non-local function Q, -P is a local func- 
tion, with P(x, £) = for |x| > 2£. In writing the above 
equations, we have dropped constant contributions to the 
tension arising from the boundaries. Indeed, if we con- 
sider a situation with periodic boundary conditions, such 
boundary terms must vanish. The integration constants 
E^ and E^P are therefore determined by the condition 
that in the final expressions l|50ll and (|5 1|) constant con- 
tributions to the tension are absent. 

The tension in the homogeneous state ^(x) = Cq, is 
given by 

X=lr,£ 3 a(c+ 2 + % 2 ) . (52) 
I 



Here, we have introduced the tension due to interactions 
between filaments of the same orientation, £={, and those 
of opposite orientation, £^±. From Eq. H46|l we obtain 

(48) 
(49) 
I 

tension is a local quantity as can be seen by introducing 
the function P(x,£) — Q(x + £) — Q{x). For the mini- 
mal model with R(x) = 1 for |x| < £/2 and R(x) = 
otherwise we obtain 



(50) 



(51) 



In this case only the interaction between parallel fila- 
ments contributes to the active part of the bundle ten- 
sion. For an oriented bundle, the tension is positive, i.e., 
contractile, whenever a > 0. 



More compact expressions can be obtained by the ap- 
proximation P(x, £) = £ for \x\ < t/2 and P(x, £) = 
elsewhere. Then 



-n j dx(J++(x) + J (x))Q(x -y) + £<J 
~j] [ dx{J+-{x) + J-+(x))Q(x -y) + 



~a?? j d£ / dx c + (x + £)c + (x)P(x -y,$)sgn(g) 



+ 2 ar] J d £ I dx c (x + £)c (x)P(x - y, £)sgn(£) 



-pr\ j d£ j dx c~ (x + £)c+ (x)P(x-y, 



Pt] \ d£ dxc + (x + £)c-(x)P(x-y,£) 



i rv+t/2 ft 

S=»(tf) w rt»7 / dx di\£,\{c+{x + i)c+(x)+c-{x + £,)c'(x)) (53) 

1 Jy-l/2 J —I 
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and 

1 ry+e/2 f £ 

Zt±(y) « -fin / di / d^(c + (a: + ^W-c1i + 0c + W) • (54) 
^ Jy-t/2 J-i 

I 



For filament distributions that vary weakly over a fila- 
ment length, this result corresponds to the expressions 
given in Ref. |22j. with fj = rj/2 and fj = ry/4. 

Contractile tension in the bundle can give rise to con- 
tractile forces exerted by the bundle. In order to il- 
lustrate this, consider a homogeneous oriented bundle 
with constant filament density c(x) = Co inside a box 
< x < L of size L while c(x) — elsewhere. In order 
to stabilize this state, we impose boundary conditions 
which immobilize filaments within the intervals [0,1] and 
[L — t, L] near the ends. Such boundary conditions could 
be realized by attaching the filaments near the end to a 
substrate. This filament distribution is stationary and 
for a < a c stable in the interval [£, L — £]. At the ends, 
the force balance Eq. JJJ is satisfied only if the force den- 
sity fext(x) = ar\£c^{l — x) is applied for x £ [0,£] and 
correspondingly at the right end. The total force acting 
on the ends is thus F = j Q f CKt (x)dx — E. This result 
indicates, that the generated force is independent of the 
bundle length and increases with the square of the fila- 
ment density. 



VI. DISCUSSION 

In this paper, we have developed a physical description 
for the dynamics and mechanics of active filament bun- 
dles. In this one-dimensional description, we describe 
the dynamics of the filament densities projected on an 
axis which is parallel to the filaments. The activity in 
these bundles results from active cross-links, that create 
relative forces between the filaments. Our approach is 
based on momentum conservation within the bundle and 
momentum exchange via external forces. Dynamic equa- 
tions can be derived most conveniently using simplifying 
assumptions such as a low motor density or low duty 
ratio, local friction, the absence of passive cross-linkers, 
and the assumption that filaments do not change their 
lengths. At low motor densities the filament currents are 
dominantly generated by interactions of filament pairs. 
Using this approach we systematically derive the minimal 
model introduced in earlier works. In general, the fila- 
ment dynamics in the active bundle is described by non- 
linear integro-differential equations. The nonlocal char- 
acter of these equations reflect the finite filament length. 

The bundle dynamics on scales much larger than 
the filament length alternatively can be described phe- 
nomenologically in a continuum limit, leading to a non- 
linear description in the form of partial differential equa- 
tions |30j. Such an approach is not limited by the above 
mentioned simplifying assumptions. However, the origin 



of the different terms which appear in the equations can- 
not be systematically related to more microscopic mech- 
anisms. Interestingly, the main features found in phe- 
nomenological descriptions are already captured by the 
more specific models described here, which are derived 
using simple approximations. This suggests that the 
main features obtained from our models can still hold in 
situations where our approximations are no longer valid. 

The equations discussed in this paper describe the av- 
erage behavior of the bundle and thus represent a mean- 
field theory where fluctuations are captured by effective 
diffusion terms but do not explicitly appear in the de- 
scription. In particular, the coefficient D of filament dif- 
fusion is effectively generated by motor-filament interac- 
tions. In the absence of such interactions, thermal fluc- 
tuations would lead to a diffusion coefficient D ~ kT/rj£, 
which becomes small for long filaments. The effective 
diffusion coefficient due to active interactions between 
parallel filaments, which do not generate a net current in 
the homogeneous state, can be estimated as D ~ Ax 2 u>. 
Here Ax is the run-length of a motor on a filament and 
oj ~ £ denotes the rate of generation of mobile cross- 
links, which grows linearly with the filament length. As 
a consequence, D ~ £ and for long filaments the diffu- 
sion is dominated by active cross- links. Numerical sim- 
ulations of computer models, which take into account 
fluctuations show, that the phenomena described in the 
previous sections persist qualitatively in the presence of 
fluctuations [22j . A thorough analysis of the effect of fluc- 
tuations will be the subject of a separate publication |30| . 

The minimal model which neglects the dynamics of the 
motor distribution already exhibits a complex scenario of 
behaviors. We discuss the full bifurcation diagram of this 
model which involves homogeneous solutions, contracted 
steady states, solitary waves and oscillatory waves. These 
states are separated by bifurcations and can partly coex- 
ist depending on parameter values. 

Taking into account the dynamics of the motor distri- 
bution does again generate the types of solutions present 
in the minimal model. Furthermore, new inhomogeneous 
steady states occur which are the one-dimensional ana- 
log of asters. In addition, the bifurcation diagrams are 
modified and are considerably richer. Contracted steady 
states are destabilized by the motor dynamics and in 
many cases become solitary waves. Furthermore, insta- 
bilities of the homogeneous state already occur in the 
absence of interactions between filaments of the same ori- 
entation, i.e. A — 0. The bifurcation diagrams exhibits 
for certain parameter ranges reentrant behavior, i.e. the 
homogeneous state is re-stabilized by increasing A or B. 

External forces modify the filament dynamics in inter- 
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esting ways. Using the momentum balance in the fila- 
ment bundle, we studied active mechanical properties of 
the bundle. In general, the bundle generates mechanical 
tension. Applying external forces to the bundle ends, a 
stationary state of a contractile bundle can be attained. 

It is interesting to compare our results to experiments, 
where purified filaments interact with motor molecules 
or small aggregates of such motors. There is one in 
vitro experiment involving filament bundles Q, where 
the contraction accompanied by polarity sorting of dis- 
ordered bundles of actin filaments in the presence of 
ATP and myosin II molecules (more specifically HMM), 
which presumably spontaneously form small aggregates, 
has been observed. Qualitatively, this corresponds to the 
state shown in Fig. [SJ However, this experiment has not 
been repeated and a systematic in vitro study explor- 
ing all the regimes discussed in this paper is lacking and 
would be very valuable. Experiments on filament bun- 
dles could also be performed using microtubules, how- 
ever, suitable preparations techniques to generate aligned 
filaments have to be developed. Artificially constructed 
kincsin aggregates would be natural candidates for mo- 
bile cross-links used in such a study. 

As mentioned above, our results are expected to be 
more general and could also apply to situations where 
additional components are present. For example it may 
be more convenient to prepare filament bundles using the 
help of passive cross-linking or bundling proteins like a- 
actinin. For a large cross-linker density or a long lifetime 
of such passive cross-links we expect qualitatively new 
behaviors. However, for low concentrations or short life- 
times of passive cross-links they are expected to mainly 
modify the effective friction coefficient and possibly other 
parameters of our model. Even in the presence of pas- 
sive cross-links, we expect the main results of our work 
to apply in certain regimes. 

A similar modification of model parameters can be ex- 
pected in general, if other proteins are present which in- 
teract with motors and/or filaments. For example, the 
interaction strength a between equally oriented filaments 
could be significantly enhanced by proteins bound to fila- 
ments which affect the speed of motors. In fact, it might 
seem odd at first glance that filaments of the same orien- 
tation exhibit significant interactions via motors at all. 
As a motor advances on both filaments with the same 
speed, no relative motion is generated. Interactions be- 
tween filaments of the same orientation are induced by 
motors that do not move with the same speed on two 
cross-linked filaments. This happens, e.g., when a motor 
arrives on one filament at the end towards which it moves. 
In this case, the motor stops on one filament while con- 
tinuing for a while to move on the second filament. This 
induces relative filament sliding via an end-effect. 

The interaction strength a between equally oriented fil- 
aments can be enhanced if the motor speed varies along 
the whole filament. For example, the speed of motors can 
be affected by the presence of other motors on the same 
filament. Such crowding would typically lead to a slowing 



down of motors as they approach the filament end |34j . 
generating relative motion of filament pairs. Further- 
more, one could imagine specific proteins bound to the 
filaments which affect the speed of motors. If such pro- 
teins had a graded distribution linked to the filament 
polarity, they would lead to strong interaction terms be- 
tween filaments of the same orientation. These examples 
illustrate, that in more complex situations where addi- 
tional proteins are present, the main results of our work 
could still apply, however, with effective parameters. 

This suggest, that the essential properties of active fil- 
ament bundles found in our work might be even more 
general and could also apply to more complex situations 
found in vivo. For example, our results could apply to 
stress fibers. These are contractile actin bundles in cells 
lacking the obvious periodic organization of muscles Q 
but containing myosins and other proteins. As we have 
demonstrated in Sect. [V] the generation of tension and 
contraction is possible through the interaction of filament 
pairs without the need of a muscle-like sarcomere struc- 
ture. The periodic boundary conditions which we use 
in several examples correspond to the situation, where 
the bundle forms a ring and could apply for example to 
contractile rings which cleaves a cell during cell division. 

Interestingly, the types of dynamic behaviors which we 
observe also include qualitatively the symmetry breaking 
presented by fragments offish keratocytes [iall3- These 
fragments consist of the lamellipodium, which is the flat- 
tened leading margin of these cells, responsible for their 
migration. Notably, they do contain neither the nucleus 
nor microtubules. These fragments exist in a symmetric 
stationary state as well as in an asymmetric locomot- 
ing state, where one can change between these states 
through sufficiently strong external perturbations |17| . 
Even though the active bundles studied in the present 
work are far from giving a description of moving kerato- 
cyte fragments, our results clearly indicate, that viewing 
the cytoskeleton as a dynamical system is a valuable con- 
cept for understanding such phenomena. Our description 
of active bundles provides a firm basis for the develop- 
ment of more profound theories of active filament sys- 
tems, which could help understanding self-organization 
and dynamic behaviors in living cells such as cell locomo- 
tion. Moving on into this direction will require a number 
of important additional ingredients. A three dimensional 
description should, for example, incorporate effects, such 
as the polymerization and depolymerization of filaments, 
non-mobile cross-linkers, capping proteins, and the inter- 
action of filaments with a membrane. 
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APPENDIX A: EIGENVALUES OF THE 
LINEARIZED TIME-EVOLUTION OPERATOR 



in Eq. (J2SJ and show, that a c > 0. 



Here, we give the complete eigenvalues of the linearized 
time-evolution operator A(fc) of the minimal model given 



The two eigenvalues of A(k) are 



J 



A ± = -k 2 + a(l - cos k)c - 2i(3kSc 

± {a 2 (l - cos k) 2 8c - (3 2 k 2 c 2 + (3 2 sin 2 k(c 2 - 8c 2 ) + 2iaf3k{l 

I 



i k)c5c} 



1/2 



(Al) 



In this expression c 



c and 5c 



The 



real part of A+, which determines the stability of the 
homogeneous state against small perturbations, is 



A(fc) = -k 2 + a(l - cos k)c + 
with 



b 2 



1 

r 



1/2 



(A2) 



APPENDIX B: NON-HOMOGENEOUS 
SOLUTIONS CLOSE TO a c 



We calculate the asymptotic solutions of the minimal 
model close to the critical value a c for arbitrary /3 and 
compare it to the expressions for small discussed in 

sect, nni 



a = a 2 (l — cos k) 2 Sc 2 

~/3 2 k 2 c 2 + (3 2 sin 2 fc(c 2 
b = 2a/3fc(l — cos k)c8c . 



6c 2 ) 



(A3) 
(A4) 



For a = this implies A(fc) = —k 2 . The derivative of A 
with respect to a is of the form 



7T = A ° 
aa 



A 1 a + A 2 a 3 



(A5) 



with Ai > 0, i = 0,1,2. It then follows that there is 
a unique critical value a c > 0, determined by A(fc = 
27r/i;a = a c ) — 0, such that the homogenous state is 
linearly stable unless a > a c and that the longest wave- 
length fitting in the system becomes unstable. 



1. Oriented bundles 



In Fourier-representation the stationary solution of 
Eq. l(2"T|) is given by 



Ok 



la ^— v cos p — 1 

k ^ p 

v 



p^k—p j 



where c_fc = c* k . In particular 



(Bl) 



J 



ci 



C'2 



2aL 2 
" 4tt 2 

ah 2 
' 4tt 2 



2tt 

cos — — 1 

ij 

2tt 

cos — 1 



cic + - 



ciCi+2 



47T 

cos — — 1 
cos — 1 

±J 



C2C 



27T 

cos — — 1 



C_iC 2 + . . 



(B2) 
(B3) 



Close to a — a c , we expand the solution in terms of c\ 
using the ansatz Cfc tx c\. Ignoring higher modes, we 
obtain 

1 ah 2 

C2 = 



1 + !^[ C0S X -l]co^ 2 
Up to third order in c\ we find 



2tt 

cos — 1 



Fia)^ - G(a)\c 1 \ 2 c 1 = 0, 



(B4) 
(B5) 



with 



F(a) = 1 



2aL 2 
4tt 2 



2tt 

cos 1 

L 



Co 



and 



(B6) 
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G(a) = 



2aL 2 
4tt 2 



47T 

cos — 1 

Li 



2tT 

cos — — 1 



The Fourier coefficients of the stationary solution are given by 

Lc 



ci = 



/ -2c 


■ 2^" 


1 — cos — 


/cosf 


Li 



1 + 1^5 [cos^-1] co- 
Ca -a c ) 1/2 + ((a- a c ) 3/2 ) 



and 



C2 



2L 2 c 2 1 1 — cos 



4tt 2 



cos 



^{a-a c ) + o({a-a c f) 



(B7) 



(B8) 



(B9) 



2. Bundles of mixed orientation 



The Fourier coefficients satisfy the following equations: 



We determine solitary waves which appear at the bi- 
furcation point a = a c using the ansatz 



c ±(x,t)=^ c ±e lfc ( 2 ^/ i+ " 



(BIO) 



2iwc 2 



live, = 



2iu>c 9 



4tt 2 . 



2a 
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C_ 1 Co 
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'2tt 



4a 



1 — cos ■ 



2tt 



2tt 



clcl + 



+4i/3 <j — c + sin — c x c{ 



1 

2 

47T 



1 — cos ■ 



47T 
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sin — c 2 c 



4tt 2 _ n 



2tt 



1 — cos ■ 



2tt 



2tt 



( c i c o - C -1 C 2 ) + 



1 — cos 



47T 
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-2i(3{ -j~c+ Cl +sin— (c+c + c\c 2 



1 4tt , _ 
-sm—cjc_ 1 + ... 



16tt 2 
I/ 2 



-c 2 + 4c? 
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1 — cos ■ 



2tt 



c i c i + 



1 — cos • 



47T 
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} 
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(Bll) 



(B12) 



(B13) 



(B14) 



We assume without loss of generality > c and where the expressions for F and G to first order in (3 are 
expand in of, which leads to given by 

F(a,/3,^)-G( a ,M| c +| 2 = 0, (B15) 



47T 2 

" L 2 



+ 2a 



2n 

1 — cos — 
Li 



cX +i 



4tt 

T /3c 



(B16) 
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and 



G = 8a- 



4tt 2 



1 — cos ■ 



2tt 



2tt 



cos ■ 



1 — cos ■ 



47T 




(B17) 



Here, A = I6n 2 /L 2 - 2a [l - cos ^] c+. 

At the bifurcation point a — a c we have F(a c , (3, uq) = 
0. In first order in a — a c we obtain for the frequency 
and the amplitude of 



uj = cj + 0Ji(a - a c ) + 



-+|2 



d a F(a c , u ) + io 1 d UJ F(a c , uj ) 
G(a Cl uj ) 



(B18) 
(a - a e ) -((BIS) 



The right hand side of the last expression needs to be 
real and positive such that we obtain the condition 



d a F(a c ,uj Q ) d u F(a c ,uj ) 



G(a c , luq) 



G(a c ,Lu ) 



(B20) 



Furthermore, the sign of the prefactor of a — a c will de- 
termine if the bifurcation is super- oder sub-critical. 
The bifurcation condition _F(a c ,o;o) = leads to 



4tt 



1 



L 2 2 [1 

47T 



cosf] c+ 



w 



-/3c, 







(B21) 
(B22) 



This implies 



d a F(a c , wo) 
9 u F(a c ,w ) 
G(a c , cj ) 



2L 2 

— [1-cos2tt£] 



'47T 2 



2tt 



(B23) 
(B24) 
(B25) 



and u>i = 0. 

The velocity of propagation is given by v — luqL/2tt = 
2/3cq~. The components c\ and c\ are given by the ex- 
pressions <|B8ll and (|B9|I for C\ and C2, respectively. For 
the distribution of minus-filaments wc find 



-2i(3 



= if3- 



sin ^c 



#-2a[l 



c 



sin ^ 1 - cos ^ 



16 
L 



27T] Cq +Cq 



c 1 



(a - a c ) + O ((t 



«c) 2 ) 



(B26) 
(B27) 



These expressions are consistent with the results of 
Sec. IIII Bl and provide analytic expressions of solitary 
wave solutions for small (3. 

APPENDIX C: FILAMENT-ADHESION 

In this section, we discuss the effects on the bundle 
dynamics due to filament adhesion to a substrate in the 
minimal model. Adhering filaments are described by the 
densities a + and a - , respectively. Attached filaments 



are assumed to be immobile, while they contribute to 
the motion of free filaments. The dynamic equations in 
dimensionless form are given by 

d t c+ = d 2 x c+ - d x J ++ - d x J + - - io a c+ + w d a\Cl) 

d t c~ = d' x c~ - d x J - d x J~ + - w a c~ + w d a~(C2) 

d t a + = uj a c + - uj d a + (C3) 

d t a~ = uj a c~ - uj d a" (C4) 

with 



J 



±± 



= d£ [ac (x 
Jo 



+ (,)+ aa ± (x + £) - ac ± (a; - £) - aer^x - £)] ^(x) 



(C5) 
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and 



d£ [j3c^( X + + 0O*(x + £)] C±(x) 



(C6) 
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FIG. 9: The effect of filament adhesion on the critical value 
etc of the minimal model. Displayed is the critical value a c 
as a function of the adhesion rate ui a for different values of 
the coupling strength /3 and for = 0.3, Cq = 0.7, a = 2a, 
[3 = 2/9, cud = 1, and L — 10. For ui a = a c = a c , which 
is the critical value in the minimal model in the absence of 
adhesion. 



Here io a and LOd are rates of attachment and detachment 
of filaments, respectively while a and characterize cou- 
pling constants between free and attached filaments. 

The homogeneous state c ± (x) = Cq and a ± (x) = 
u a c /uid is stationary. It becomes unstable at a criti- 
cal value a c which for oriented bundles is given by 



' fJ,LV a + UJ d ' 



(C7) 



where a c is the critical value in the minimal model with 
uj a = and for = a/a. For ^ a Hopf- bifurcation 
occurs, the dependence of the critical value d c on the 
attachment rate uj a is shown in Fig. [5] for a = 2a and 
= 20. 

In the case — 0, the stationary distributions are 
given by (^{x\u) a ) = j^f^c^ (x; 0) and a ± (x;uj a ) = 



uu^+uT c ( x 'i Q) wnere c ± (x;0) denote the stationary 
states of the minimal model. For 0^0 solutions can 
be obtained in the limit of small attachment rates by 
expanding around solitary waves for uj a = 0. Consider a 



solitary wave in the minimal model given by c ± (x, t; 0) = 
u o ( x ~ v ot) aud a ± (x, t; 0) = 0. Assuming that for finite 
uj a solutions of the form c ± (x, t; ui a ) = u ± (x — vt) and 
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FIG. 10: The total filament current I of propagating filament 
patterns with filament adhesion as a function of the adhesion 
rate u a . Symbols represent numerically obtained results for 
different detachment rates u>d- Parameters are a = 1.5, /3 = 2, 
$ = 2/3, cj = 0.3, cjT = 0.7, and L = 10. The maximal 
current decreases and is reached at larger ratios of uj a /ujd as 
u>d is increased. 



a ± (x,t : Lu a ) = r (x — vt) exist, we write 



v = v + uj a vi + 

U ± = Uq + LV a uf -| 

= War^ + . . . 



(C8) 
(C9) 
(C10) 



The equations for can be solved explicitly. We find 



- I dt' e Udt 'u±{x-v t') , (Cll) 

1 ./n 



where T = vqL. 

This result can in turn be used to calculate in first 
order in co a the total net current I associated with these 
solitary waves. Indeed, in lowest order this current is 
given by 
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I = 2Lo a \ dx\aJ^d([(4(x + Z)-rf(x-0)ut(x) + (ri(x + 0-rT(x-£))uo(x)] 

+/? [ d£ [rt(x + Ouo(x)-ri(x + OuUx)]} > (C12) 



where, for simplicity, we have chosen a = 2a and f3 = 2/3. 
Evaluating this expression using Eq. (|C11|) and u J (z) = 



Er=-co4«e i2 ™ /L we obtain 



I = 8^ a L 3 Uo a]T A n (\u+J 2 + | u -„| 2 ) + 0(ul) , 

n=l 

(C13) 



where 



cos ^ - 1 
2 L 2 + 47r 2 n 2 w 2 



(C14) 



r 



with wo and depending on j3 as described in Sec. lIIIBl 
This result shows, that solitary waves are accompanied 
with a net filament transport if filament adhesion occurs. 
Since A n < for all n, this transport occurs into the 
opposite direction as wave propagation. Remarkably, as 
indicated by the prefactor a, it is the interaction between 
parallel filaments that generates the current. Note, that 
for (3 = we have vq = and thus 1 = 0. Numerical 
results for / as a function of u a are shown in Fig. 1101 
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